clear
set mem 80m
set more off
set matsize 800
capture log close

global citylist "d_Chengdu d_Chongqing d_Haerbin d_Hefei"

cd "D:\MyPapers\Happiness\TWINSDATA\code\2023May-Dataverse"
use nontwin_PNAS.dta, clear 

//Table 1 (summary statistics), Column 4: NBS Non-twin sample
sum happy_4s income_num age male  birth_weight VeryEarly_disease_d early_disease edu working work_over_40h married divorced widowe population 

//Table S4: OLS with NBS Non-twin sample
quiet reg happy_4s income age agesquared male $citylist,robust
outreg2 using TableS4_OLS_NonTwins.doc,  nolabel replace coefastr se drop(d_*)

quiet reg happy_4s income age agesquared male  birth_weight VeryEarly_disease_d $citylist,robust
outreg2 using TableS4_OLS_NonTwins.doc,  nolabel append coefastr se drop(d_*)

quiet reg happy_4s income age agesquared male  birth_weight early_disease_d edu $citylist,robust
outreg2 using TableS4_OLS_NonTwins.doc,  nolabel append coefastr se drop(d_*)

quiet reg happy_4s income age agesquared male  birth_weight early_disease_d edu working work_over_40h $citylist,robust
outreg2 using TableS4_OLS_NonTwins.doc,  nolabel append coefastr se drop(d_*)

quiet reg happy_4s income age agesquared male  birth_weight early_disease_d edu working work_over_40h married divorced widowed $citylist, robust
outreg2 using TableS4_OLS_NonTwins.doc,  nolabel append coefastr se drop(d_*)

quiet reg happy_4s income age agesquared male  birth_weight early_disease_d edu working work_over_40h married divorced widowed population $citylist,robust
outreg2 using TableS4_OLS_NonTwins.doc,  nolabel append coefastr se drop(d_*)

